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A cold electron beam propagating through a background plasma is subject to filamentation process 
due to the Weibel instability. If the initial beam radius is large compared with the electron skin depth 
and the beam density is much smaller than the background plasma density, multiple filaments merge 
many times. Because of this non-adiabatic process, the beam perpendicular energy of initially cold 
beam grows until all filaments coalesce into one pinched beam with the beam radius much smaller 
than initial radius and smaller than the electron skin depth. It was shown through particle-in-cell 
simulations that a significant fraction of the beam is not pinched by the magnetic forces of the 
pinched beam and fills most of the plasma region. The resulting electron beam energy distribution 
in the perpendicular direction is close to a Maxwellian for the bulk electrons. However, there are 
significant departures from a Maxwellian for low and high perpendicular energy (deeply trapped 
and untrapped electrons). An analytical model is developed describing the density profile of the 
resulting pinched beam and large low-density halo around it. Based on this analytical model, a 
calculation of the energy transfer from the beam longitudinal kinetic energy to the transverse beam 
kinetic energy, the self-magnetic field, and the plasma electrons is performed. Results of analytical 
theory agree well with the particle-in-cell simulations results. 



PACS numbers: 52.35.-g, 52.35.Py, 52.35.Qz 



I. INTRODUCTION 



The Weibel instability (WI) [1-5] of plasmas with 
anisotropic velocity distribution is one of the most basic 
and long-studied collective plasma processes. For exam- 
ple, propagation of an electron beam through the back- 
ground plasma is subject to strong WI [2|. There has 
been a significant revival in theoretical studies of the WI 
because it is viewed as highly relevant to at least two 
area of science: astrophysics of gamma-ray bursts and 
their afterglows and the Fast Ignitor [l2j scenario 

of the inertial confinement fusion (ICF). Specifically, gen- 
eration of the upstream magnetic field during the GRB 
aftershocks is considered necessary for explaining emis- 
sion spectra of the afterglows as well as for generating and 
sustaining collisionless shocks responsible for particle ac- 
celeration during GRBs. Collisionless Weibel Instability 
is the likeliest mechanism [rjl-fll| for producing such mag- 
netic fields. The WI is likel y to play an important role 
in the Fast Ignitor scenario [l2[ because it can result in 
the collective energy loss of a relativistic electron beam 
in both coronal and core plasma regions |1 2 - 20]. Be- 
cause the relativistic electron beam has to travel through 
an enormous density gradient (varying from 10 22 cm -3 
near the critical surface where the beam is produced to 
10 26 cm~ 3 in the dense core), both collisionless and colli- 
sional WI manifest themselves along the beam's path. 

The dynamics and energetics of the nonlinear satura- 
tion and long-term behavior of the WI are important for 
both laboratory and astrophysical plasmas. For example, 
collisionless shock dynamics depends on the long-term 



evolution of the magnetic field energy. Specifically, it is 
not clear whether the long-term magnetic fields gener- 
ated during the coalescence of current filaments remain 
finite fltl ] or decay with time [|| (and if they do, according 
to what physical mechanism). Numerous numerical sim- 
ulations jjj [2l| demonstrated that magnetic field energy 
grows during the earlier stages of the WI and starts de- 
caying during the later (strongly nonlinear) stage. The 
reason for this decay has never been fully understood. 
One decay mechanism based on the merger of filaments 
bearing super- Alfvenic current {I > I a = jf3mc 3 /e [22J, 
where — e and m are the electron charge and mass, re- 
spectively, and c is the speed of light in vacuum) during 
the late stage of the WI has been recently identified [23| . 
In the review [24| we presented a detailed analysis of the 
high-current filaments' current and density profiles and 
provide qualitative and quantitative explanation of the 
energetics of their merger in the limit (7 > I a). In that 
case the filaments carry super Alfenic current and their 
perpendicular energy distribution is closely described by 
a RH distribution that is 5— function of perpendicular 
kinetic energy, or KV distribution as it is called in ac- 
celerator physics. For such a distribution functions the 
density profile in the filament is flat and the radial electric 
field vanishes. Particle-in-cell simulation of the nonlinear 
stages of the Weibel instability showed significant ion ac- 
celeration in the radial electric field. Therefore, it is 
very important to investigate departure from RH 
or KV distribution during nonlinear stages of the 
Weibel instability, which ultimately determines in 
acceleration in the micro field of the filammetns. 

To that end we investigated propagation of the rela- 
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tivistic electron beam through a background plasma and 
development of the self-electric and magnetic field. Of 
particular interest is a case when the beam transverse size 
is much larger than the electron skin depth by a factor of 
ten and more, r b > I05 p . The filament resulting from the 
Weibel instability are typically of the si ze of the ele ctron 
skin depth, S p . Here, S p = c/u> p , u> p — y/ 'Ai:e 2 np/m is the 
electron plasma frequency, and n p is the uniform back- 
ground plasma density. Therefore in the limit r b > 106 p 
many, > 100 filaments are formed and then merge mul- 
tiple times. In the process of merger there is an effec- 
tive exchange in the perpendicular energy between beam 
particles and plasma electrons. If the beam density ap- 
proaches the background plasma density, the beam elec- 
trons expel the plasma electrons and the beam electron 
space charge is neutralized by the background ions. In 
this case the further pinching of the beam is limited by 
the background plasma density. To avoid this limitation, 
we studied the very low density beam with the density a 
factor of 1000 less than the background plasma density. 
This insured that as beam filaments merge and the beam 
density dramatically increases the beam density still re- 
mains small compared with the plasma density. In this 
limiting case, common PIC codes are not efficient for the 
description of the plasma electron due to the large nu- 
merical noise compared to the beam density. However, 
we can utilize a semi-analytic approach for description 
of the plasma electrons. Because the beam evolution 
occur on a time scale much smaller than the electron 
plasma period, the electron background plasma adiabat- 
ically modifies to the beam current profile via the return 
current, see Ref. [24| for details. We also assume charge 
neutrality of the system consisting of the beam electrons, 
ambient plasma electrons, and ambient plasma ions fl9| . 
These simulations do not resolve the motion of the am- 
bient plasma electrons and, therefore, can take compu- 
tational time steps At > l/ui p . After exclusion of fast 
motions of the plasma electrons, the beam particles are 
treated as macroparticles in PIC algorithm. 



A. Low-noise efficient quasi-neutral particle-in-cell 

code 

The logic behind the quasi-neutral code is that the full 
dynamics of the ambient plasma need not be simulated, 
and its density can be obtained from the quasi-neutrality 
condition: 

n pe (x) = n - n b {x). (1) 

Therefore, ambient plasma is modelled as a passive fluid 
that responds to the evolving electron beam in order to 
maintain charge neutrality. Electron beam particles are 
modelled using numerical macro-particles that are ad- 
vanced in time by the self-consistently determined elec- 
tric and magnetic fields. The leading magnetic field 
B±_ = —e z x Vj_"0 develops in the x — y plane, where 



if) is the z-component of the vector potential. The induc- 
tive electric field associated with the time- varying flux ip 
is E z — —(l/c)dtip. Electric field also has a transverse 
component E±_ that is found from the quasi-static force 
balance of the ambient plasma electrons in the plane: 
E±_ + v pz x B±/c — 0, where v pz = v pz e z is the return 
flow of the ambient plasma. This quasi-equilibrium is 
the consequence of another observation from direct PIC 
simulations: that the transverse velocity v p ± of the am- 
bient plasma electrons is considerably smaller than the 
beam's average transverse speed \v b ±\ and plasma's lon- 
gitudinal velocity v pz . One of the consequences of that 
is that the dominant magnetic field is the transverse 
one [33j], i. e. that the out-of-plane magnetic field is small: 
|V"0| » B z . To summarize, these are the dominant elec- 
tric and magnetic fields of the quasi-neutral beam-plasma 
system: 

- 1 dip - 

B± = -e z x V±ip,E z = ---Z l E± = -{v pz /c)V ±_ijj. 

(2) 

For collisionless plasma, two important conservation 
laws simplify the description of the plasma motion: con- 
servation of the canonical momentum in the z-direction 
and the conservation of the generalized vorticity [l9|, H3, 
[3B|. The former is essential for deriving the field equa- 
tion for ip that defines the dominant in-plane magnetic 
field. Conservation of the canonical momentum trans- 
lates into the non-relativistic expression for the plasma 

return velocity: v pz /c = ip (or v pz /c — ip/ \J \ + ip 2 in the 

rclativistic case), where ip — eip/mc 2 is the dimensionless 
vector potential. From Ampere's law then follows that 

V\ip - 4irenv pz /c — -4irJ bz /c, (3) 

where the displacement current is neglected to be con- 
sistent with the quasi-neutrality assumption. The z- 
componcnt of the beam electron current is calculated 
from the beam's macro-particles' contribution: J bz — 
— qjVj Z 5 2 (x '— Xj) I 'L z , where index j labels numerical 
macroparticles, and Qj/L z and Mj/L z are the macro- 
particle's charge and mass per unit length. For non- 
relativistic collisionless plasma electrons Eq. ([3]) is sim- 
plified to Vj_t/> — k pe tp = —4irJi, z /c, where the spatially- 
nonuniform k pe = u pe /c 2 is obtained from Eq. (fT]) 

through oj pe (x) — iire 2 n e (x)/m. Specifically, Vj_ X Jb± 
is small in the linear limit n b << n p 

As in the standard PIC, beam electrons are mod- 
elled kinetically using macro-particles with the effective 
per-unit-lcngth charges and masses qj and rrij satisfying 
qj/rrij — e/m, where index j labels numerical macro- 
particles. The longitudinal momentum of a beam elec- 
tron (assumed collisionless owing to its relativistic en- 
ergy) is found from the conservation of the canonical mo- 
mentum: 

li v i* = IjoVjto + — (tp - ipjo), (4) 
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FIG. 1: The beam density rib l n w in the beam cross-section at 
different moments of time. First, the beam breaks into a large 
number of non-stationary filaments in which fluctuations of 
the beam density are of the order of ribo- Non-stationary fil- 
aments eventually coalesce into one quasistationary filament. 



where we assume that the initial field in the plasma van- 
ishes: ipjo = 0. The transverse equation of motion for 
the beam electrons is: 



d('jjVj±) _ e(Vj 



dt 



(5) 



where the second term in the rhs of Eq. ([5]) is due to 
the extra pinching of the electron beam provided by the 
transverse ambipolar electric field E±_ = — (v pz / c)V ±ip 
that develops in response to the v x B expulsion of the 
electron plasma fluid. Note that Ej_ counters the mag- 
netic expulsion of the ambient plasma, yet reinforces the 
magnetic pinching of the beam. The above expression for 
E± are only valid in the absence of the complete plasma 
expulsion from the beam filaments, and needs to be mod- 
ified when such expulsion takes place. 



II. STAGES OF THE WEIBEL INSTABILITY 

During the extensively studied [j| [H, HH early 
stage of Weibel instability, the electron beam of the den- 
sity ribo < rip and the radius i?; n it 3> c/uj pe breaks up 
into a large number of filaments, see Fig. [TJ We model 
low-density electron beams with large cross-section ra- 
dius by performing simulations in square-shape domain 
using periodic boundary conditions (in the simulations, 
we chose n,bo/n p = 0.001 and k p L = 10). Small random 
perturbations are imposed initially on the homogeneous 
beam density. The beam temperature is assumed to be 
small, the relativistic factor 7^0 = 10, and the beam cur- 
rent is compensated by the plasma return current. 
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FIG. 2: Evolution of the magnetic energy spectrum of per- 
turbations. Growth of short-wave perturbations is suppressed 
due while long-wave perturbation grow until the system 
reaches equilibrium. 




FIG. 3: The evolution of the transverse kinetic energy of beam 
particles (blue), the magnetic energy (green), and the energy 
of plasma particles (red). The saturation level are of Non- 
zero initial value of the plasma energy is determined by the 
energy in the return current. The energy is measured in units 
NT,, where TV is the number particles in the beam and T, = 
mc^nbo/rip. 



When perturbations are small, the magnetic field energy 
grows exponentially with the growth rate of order 7wi 
given by 



7wi(fc±) = u p 



kp 



(6) 



where f3 x — v bz /c and k± is the transverse component of 
the wavevector. Although all modes with k±_ » k p grow 
simultaneously (with the rate 7w = Jwi(k± —> 00)), 
wavebreaking occurs first in short-wave modes and then 
in long-wave modes. Subsequent formation of multi- 
stream flow in the transverse beam motion effectively 
increases the transverse temperature of the beam par- 
ticles. The thermal motion, in turn, suppresses growth 
of small-scale perturbations while large scale perturba- 
tions continue to grow. The dynamics of this complex 
process is illustrated in Fig. [2] where the evolution of 
the magnetic energy spectrum of perturbations is pre- 
sented. The evolution of the transverse kinetic energy 
of beam particles and the magnetic energy is given in 
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FIG. 4: The beam particle distribution ln[/(ex)//(0)] over 
transverse kinetic energy e± (measured in units T»). The 
transverse temperature extracted from this dependence T w 
1.2T». 



Fig. [3J In the course of multiple merging, non-stationary 
filaments become bigger and the distance between them 
increases. Eventually, they coalesce in one quasistation- 
ary filament (in our example, it happens at t ~ 
Often, this filament has elliptical form and slow rota- 
tion is observed. Particles in the filament perform com- 
plicated oscillatory motion along nonclosed trajectories. 
Nevertheless, our simulations show that the distribution 
function of the trapped particles over the transverse ki- 
netic energy is approximately Maxwellian, see Fig. [4] Ap- 
parently, Maxwellization is caused by multi-generation 
merging of filaments. 

Note that due to wavebreaking, fluctuations of the 
beam density become significant (of the order of ribo) 
almost from the very beginning. As seen from Fig. [TJ 
the density in the quasistationary filament exceeds tens 
times the initial level of the density. Although it is not 
noticeable in the figure, there is a big fraction of beam 
particles (in our case about 30%) outside of the filament 
in the large area where the final beam density is small, 
Tib < n-bo- These particles get dispersed in the process of 
multiple filament merging, and gain sufficient transverse 
temperature forming stable (against further development 
of Weibel instability) background . 



III. ESTIMATES OF THE QUASISTATIONARY 
FILAMENT PARAMETERS. 

Characteristics of the quasistationary filament can be 
estimated by the following way. Making a rough as- 
sumption that all particles from L x L area are gath- 
ered in small area of the radius of the order of skin- 
depth length S p = kp 1 , one can estimate the beam 
density there as rib ~ ^fco * k 2 L 2 and the magnetic 
field as B ~ Anerib/kp, so that B 2 ~ 47rmc 2 rc. 2 /n p ~ 
Anrnc 2 {n 2 Q /np)k 2 L 2 . Equating the transverse thermal 
pressure in the filament to the pressure of the magnetic 
field, nbT± ~ B /8ir, one can estimate the transverse 



temperature in the beam filament 

T± ~ mc 2 ^k 2 v L 2 ~ 2irmc 2 -^, (7) 
2n p F Iao 

where Iao = ftmc 3 /e. This is a very rough estimate. 
More accurate formula (|27p which takes into account 
the spatial distribution of beam particles in the fila- 
ment gives considerably smaller numerical coefficient. 
The fraction of the longitudinal kinetic beam energy, 
W\\ = 7f,mc 2 (nboi 2 ), transferred to the transverse beam 
motion is given by 

'^pI±, (kr Lf(^) 2 , (8) 
W\\ \uj p ) 

The energy transferred to the magnetic field and to 
the plasma is of the same order. Since the beam den- 
sity rib in our consideration is always less than ri p , fil- 
aments of the radius c/lo p carry sub-Alfvenic currents, 
enf,c/3 z 7rc 2 /Wp << f3 z mc 3 /e. 

IV. SIMILARITY AND CONSERVATION LAWS 
FOR THE BEAM FILAMENTATION DYNAMICS 
IN THE LIMIT OF LOW DENSITY 

When rib << %> during entire beam evolution and 

Vb ~ (0,0, —c), for analytical tractability we introduce 
further simplification in our hybrid model. To the first 
order in the small parameter nb/n pi the beam dynamics 
is governed by the equations: 

k~ 2 Aip = ip - (n b /np)mc 2 /e, (9) 

From these equations one can find that the evolution of 
the low density beam beam obeys a similarity law: it re- 
mains the same when space coordinates, time, the beam 



density and vector potential rescale as 

rj. oc A p , (11) 

t^ 1 cx ujp^nbo/^brip — jw, (12) 

v± oc cy/ribo/jbrip — c(7vr/w P ), (13) 

n b (x n b0 , (14) 

ip cx (n b0 /n p )mc 2 /e, (15) 

7t,T ?-^ oc mc 2 nb/rip. (16) 



In particular, it means that at small beam densi- 
ties it is enough to study beam propagation only 
at the one level of the initial density. The propa- 
gation at other levels can be found by rescaling the sim- 
ulation results according to Eqs. (fTTj) - (fl"5j) . As seen from 
Eqs. jH]), (US]), (151) and (IB]), big beam relativistic fac- 
tors 7fc slow down the beam filamentation not changing 
magnetic fields created by filaments and transverse ki- 
netic energy of beam particles. 
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It follows from the conservation of longitudinal 
momentum that for ultrarclativistic beam 7&mc 2 ~ 
7bo'7ic 2 + m'*fbo'mv\/2 — e(tp — -0*), where ip* is the initial 
value of the magnetic potential. The energy conservation 
law 13 711 can be now transform to the form 



If 1 

W± = - n b mj b v 2 L dS - - 



eribipdS = const. (17) 



This formula can be also derived directly from Eqs. Q 
and (fTU|). The potential energy is given by the second 
term in the right-hand side of Eq. ([TT]) : it includes the 
energy of the magnetic field and the plasma motion: 




5 
k y 



FIG. 5: Crosssections of the filament density (blue) in x 



and y directions at the moment t w 5A-y v 



Green lines 



correspond to Bennett pinch with the radius tb = 0Ak~ . 
The temperature of the filament is taken from simulations, 
T — l.35mc 2 ribo/np- 



eribipdS = — I B 2 dS 



1 

8^ 



n T1 mv 2 dS, 



(18) 



where B = V?/> and v — eip/mc. The coefficient 1/2 
in front of the potential term reflects the fact that each 
particle is counted twice in the integral (|T7|) so that the 
total energy of the one particle (which in general varies 
with time!) is 



eip. 



(19) 



The magnetic potential forms a potential wells for parti- 
cles, U = —eip; these wells merge with each other when 
corresponding filaments merge. 

The other obvious integral on the motion of the system 
([9])- ([TO]) is the conservation of the total number of the 
macroparticles: 



N 



ribdS = const. 



(20) 



Note that for the beam with initially negligible transverse 
velocity spread homogeneously distributed over the area 
L x L, the magnetic potential ip — (rib / n p )mc 2 / e and 
the total energy and number of particles are given by 
formulas, 



W± 



1 



-Neip* 



1 



-NT*, N = n b0 L 2 , 



(21) 



where eip* = T* = (ribo/n p )mc 2 . 

It is instructive to note that, in the framework of 
Eqs. and (fT0|) , the magnetic potential can be pre- 
sented as a sum of potentials created by each beam par- 
ticle, 



iP(r,t) = 



2 , nen ri 



K (k p \r-r>\)n b (f*,t)df*. (22) 



Thus, the interaction between particles in the low density 
beam is a pairwise Coulomb attraction screened by Bcsscl 
function at distances larger than the skin-depth A" 1 , and 
the potential energy of beam particles can be presented 
as a sum of potential energies of all particle pairs. 



V. STRUCTURE OF FILAMENTS 

Figure [5] shows vertical and horizontal density cross- 
sections of the beam corresponding to the density distri- 
bution in the quasistationary filament in Fig. [T] As can 
be seen from this figure the filament pinched to the radius 
even smaller than the skin depth S p = k p l . For such ra- 
dius the plasma screening of the magnetic field does not 
occur. Therefore, the magnetic pinching force is balanced 
by the pressure gradient Vp — j x B and V x B = Airj/ c. 
This equilibrium correspond to the Bennett pinch |22j 



n b 



S/(k p r B ) 



T 



(1 + (r/rg) 2 ) 2 mc 2 



(23) 



For the Bennett pinch the self-magnetic field outside 
of the filament is decreasing with radius B = Airl/cr and 
magnetic potential ip decreases as — 4tt(I/c) lnr. There- 
fore according to the Boltzmann relationship the density 
n b * exp(eip/T) decreases as power law of radius. Due to 
plasma screening of the magnetic field, the magnetic field 
and variation of the magnetic flux vanishes at r » S p . 
Therefore density does not approach zero at large radius 
but tends to a finite value. This means that a low den- 
sity halo forms outside filament. Detailed analysis shows 
that the total number density of particles in the halo is 
comparable to the total number density in the filament. 
See detail calculation in Appendix. Difference between 
Bennett distribution and modified Bennett distribution 
is illustrated in Fig. O 




FIG. 6: The beam particle density in Bennett pinch (red) and 
modified Bennett pinch (blue), k v rs = 0.6. 

The energy distribution of beam particles in the phase 
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correspond to the straight line). However, there is sig- 
nificant departures from a Maxwellian for low and high 
total energy (deeply trapped and untrapped electrons). 
This because these electrons do not undergo many en- 
ergy transfer with other electrons during filament merger. 
Low energy electrons are confined to the center of the fila- 
ment and do not experience the time dependent magnetic 
field and inductive electric field during filament merger. 
Similarly untrapped electrons do not participate in en- 
ergy transfer. Therefore, tail of the EEDF is strongly 
depleted at e tot > i/Wx - tp m m ■ 



FIG. 7: Contour lines of the beam particle distribution 
ln[///(0)] in the space formed by the transverse kinetic en- 
ergy e± = 7i,m«x/2 and potential energy U = —eip + eip max 
(measured in units T* = mc? / n p ) . The variable U changes 
from zero (in the center of the filament, where eip — eipmax ~ 
7.6mc 2 n(,o/iip) to the value 7 .3mc 2 'n^o / 'n p outside of the fil- 
ament (where eip = eip m i n sss 0.25mc 2 nbo/n p ). Separatrix 
(black line) divides the phase space into two regions: the 
first region (in the shape of triangle) corresponds to parti- 
cles trapped in the filament and the second one to untrapped 
particles. The point e± — 0, U = corresponds to the cold 
particles at the center of the filament. 



1-2.5 
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FIG. 8: The beam particle distribution ln[///(0)] as a func- 
tion of etot (measured in units T*). Dashed line corresponds 
to the distribution function when energy exchange between 
electrons is absent. 



space space (e _l , U) obtained from our simulations is 
shown in Fig. [7l (For sake of convenience, the poten- 
tial energy here is defined as U — —eip + eip max .) One 
can see from this figure that the distribution function 
depends approximately only on the total particle energy 
£tot — + U. Such dependence is a simple consequence 
of the phase mixing of particles moving in the quasista- 
tionary potential well along trajectories with close total 
energies. 

The electron distribution as a function of etot is shown 
in Fig. [5] Note that if there is no energy exchanged be- 
tween electrons, the beam particle distribution / over 
total energy e tot trapped in a potential well is constant 
[35l ] . Because of the many events leading to the en- 
ergy exchange during filament merger the energy elec- 
tron distribution function is close to a Maxwellian (for 
ideal Maxwellian distribution function, ln[///(0)] should 



VI. ENERGY TRANSFER FROM THE BEAM 
LONGITUDINAL KINETIC ENERGY TO THE 
BEAM TRANSVERSE KINETIC ENERGY SELF 
MAGNETIC FIELD AND PLASMA ELECTRONS 

From two conservation laws, conservation of the num- 
ber of particles and the energy conservation law, one can 
find the parameters T and of the modified Bennett 
pinch and then the energy transfer from the beam lon- 
gitudinal kinetic energy to the beam transverse kinetic 
energy, magnetic field and plasma electrons can be cal- 
culated. 

Using the Boltzmann distribution of the density in qua- 
sistationary state, the conservation of the number of par- 
ticles gives 



n b * Jexp(etp/T) dS = N, 



(24) 



where TV = n^L 2 . For Maxwellian distribution function, 
the energy conservation law (1171) with the initial constant 
(f2"T|) can be transformed to 



TN = -n u 



exp(eip/T)e(ip - ip*)dS, (25) 



where ip* is the initial value of the magnetic potential. 
The term in the left-hand side corresponds to the trans- 
verse kinetic energy of the beam particles, and the term 
in the right-hand side corresponds to the increase in the 
effective potential energy. This increase is the same as 
the increase in the energy of the transverse motion (that 
is, < Sip >— T). Thus, the energy drawn from the longi- 
tudinal motion of the beam during its evolution is: 



(26) 



< A(7 b TOc 2 )|, >« -2T, 



After substituting the magnetic potential from Eq. (|47)) 
into Eqs. (|24p and (|25p and performing integrations, one 
can find the following asymptotic formula 



T 



k 2 L 2 



T, 



mc 2 /2 



1 



W{a klL 2 ) 8tt 



W(a k 2 p L 2 ) I A0 ' 



(27) 



where W(x) is the product logarithm (or Lambert W- 
function defined by the equation VFexp(VF) = x), a = 
exp(2 7i5 - 2)/4tt = 0.0341... and 7^ = 0.5772... is 
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the Euler constant, T» = mc 2 n b o/n p and h/lAo = 
en bQ cL 2 /(mc 3 /e). For fc p L = 10, Eq. (gTJ) gives T « 
1.7T* while simulations show that the transverse elec- 
tron energy averaged over all beam particles is 30% less: 
< e± >= T ps 1.2T*. The temperature is higher in the 
filament T sw 1.4T* and lower outside of it. Difference 
between theoretical value of T and that obtained from 
simulations is due to incomplete Maxwellization of the 
beam electron distribution function. 

Since the product logarithm W(x) ~ \n(x/lnx) , one 
can conclude from Eq. (|2~T1) that the transverse temper- 
ature of the low-density beam slowly decreases as the 
beam cros-ssection aria increases (while the beam cur- 
rent is kept constant, = const). 

At large beam cross-section, the energy transfcred from 
the beam longitudinal motion to the beam transverse mo- 
tion, the magnetic field and the plasma motion is divided 
in the following way: 



n b0 <e±> _ 1 
2n b0 T 2 > 



<B 2 /&rr> ^ 
2n b0 T ~ 

<Ti p Ae p > ^ l 
2n b0 T 2 



-In r B + 1.3rfl-0.88 
l + W(a klL 2 ) ' 

- lnrB+1.3rB-0.88 
l+W{a klL^) ' 



(28) 
(29) 
(30) 



For k p L = 10, Eq. fl52J gives r B w 0.4 and Eq. (|29[) 
gives < B 2 /8ir >/2n bQ T w 0.26, and Eq. $U$) gives 
< n p Ae p >/2n b0 T « 0.24. One can find from Fig. [3] 
that the saturation level of the transverse kinetic energy 
of the beam (1.2NT*) is two times higher than the sat- 
uration level of the magnetic energy (0.6 AT*). It is also 
two times higher than the saturation level of the plasma 
energy without initial energy (0.6 AT*). This observation 
is consistent with theoretical prediction. 



VII. CONCLUSION 

We have considered the filamentation process of an 
electron beam with small density and large radius prop- 
agating through the dense plasma. We have found that 
the beam dynamics obeys simple similarity laws which 
simplify the beam description. During the formation of 
filaments and their merging the total transverse energy 
of the transverse motion of the beam particles is con- 
served. The effective transverse potential energy is com- 
prised mostly of the energy of the magnetic field and the 
plasma motion and can be presented as a sum of pairwise 
interaction energies between beam particles. 

Our particle-in-cell simulations have shown that all 
filaments eventually coalesce into one pinched beam al- 
though a significant fraction of the particles remains un- 
trapped and scattered outside of this pinch. The electron 
beam distribution over the transverse kinetic energy is 
close to the Maxwellian one for the bulk of electrons. We 
have developed an analytical model and found the distri- 
bution of the particles in the modified Bennett pinch and 
in the low-density halo around it. In particular, we have 
found that the radius of the Bennett pinch modified by 



the return plasma current is always less or of the order of 
the plasma skin depth. For a given beam transverse size, 
we have calculated the energy transfer from the beam 
longitudinal motion to the transverse kinetic energy of 
the beam, the self-magnetic field and the plasma motion. 
It has turned out that the final transverse beam temper- 
ature is proportional to the beam current and slowly de- 
creases with the beam radius. The results obtained from 
the model agree well with those obtained from particle- 
in-cell simulations. 

Note that although our consideration is limited to ini- 
tially cold beams, the calculations can be performed in a 
straightforward manner for the arbitrary case. 



VIII. APPENDIX I. DETAILED 

CHARACTERISTICS OF MODIFIED BENNETT 
PINCH 

In this section we describe the structure of an iso- 
lated filament similar to one shown in Fig. [TJ We as- 
sume that all particles have Maxwellian distribution over 
velocities, i.e. the system is in a thermodynamic equi- 
librium. Hence, the density of particles is expected to 
be distributed in the effective potential tp according to 
Boltzmann law: 



rib = n b * exp(e^/T), 



(31) 



where n b * is some constant and T is the transverse tem- 
perature (equal to the average transverse energy of a 
beam particle T —< m^ b v\j2 >). 

Substitution of Eq. (l3~Tj) into Eq. © gives closed equa- 
tion for the vector potential 



— n 6 *exp^— j 



(32) 



Since the magnetic filed is small far from the filament, 
the magnetic potential should be constant at large dis- 
tances from the filament. Under this condition, equation 
(|32f determines the function ip{x,y) at each value of the 
parameter n b * . In one dimensional case when all quanti- 
ties depend on one transverse coordinate x, one can solve 
Eq. (|3"2"j) analytically and find the potential and density 
distribution in an isolated filament, see Appendix II. 

In two dimensional case when particles are gathered 
from large area LxL ( much bigger than the aria occupied 
by the filament), the fully thermalized filament is axially 
symmetric and therefore Eq. (|32[) can be transform to 



1 d d /eijj 

l^d^ r d^\Y 



eip 



(f-^)-(33) 



where ru = [T /A-Ke 2 n b (0)] 1 ^ 2 is the Debye radius cor- 
responding to the beam density at the filament center, 
n b (0) = n b * exp(eipo/T) 1 and ipo = ^(0)- This equation 
can be solved numerically by the shooting method: at 
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given value of the parameter k p rr>, we find such eipo/T 
at which the magnetic potential tp is a monotonically 
decreasing function of r with asymptote ip — > ipoo and 
dip/dr — > at r — > oo. 

It turns out that the solution with these properties 
exists only when k p ro < 1. Indeed, at the point r = 
where the magnetic potential reaches maximum, we have 
%p = ipo, dip/dr = 0, and d 2 ip/dr 2 < 0. Therefore, the 
right-hand side of Eq. (|33|) should be negative at this 
point, that is, 




eipo 
T 



< 



(34) 



On the other hand, the expression in the right-hand side 
of Eq. (|33|) should be equal to zero at large r. Therefore, 



V 



FIG. 9: The magnetic potential at different values of the 
Debye radius: k p ro — 1 (black), k p ro = 0.23 (red), and 
k p r D = 0.022 (blue). 



eip c 



T 



— j- exp [~7^~ jrj > ~jr ex P [~7^~ — "^^35ljhis equation describes the Bennett distribution of beam 

'•' " ' ' ' particles with uniform transverse beam temperature 22] : 



p D 



After simple manipulations, we find that i/'oo and ipo sat- 
isfy to the following condition, 



F 



eip Q 



T 



< F 



T 



ei/'oo ^ eipo 
at ~y~ < — > ( 36 ) 



where F(x) = — a;exp(— x). Since this function grows 
only at x > 1, F'(x) = (x — l)exp(— x), the condition 
(|36[) imposes restriction on the maximum value of the 
magnetic potential, etpa/T > 1. It means that filament- 
like solutions exist only when l/k 2 r 2 D > eipo/T > 1, that 
is, only when k p ro < 1. Note that the trivial solution 
of Eq. (l33|) eip/T = const = l/k 2 r 2 D exists at all values 
of k p ro- However, when k p rr> > 1, there are no other 
solutions except the trivial one: the Weibel instability is 
suppressed by the thermal motion of particles, and the 
beam remains homogeneous with n b = const. 

Figure [5] shows the magnetic potential as a function 
of r at different Debye radii. Note that eipo/T = 1 at 
kpro — 1, eipq/T « 4.3 at k p r£> = 0.23, and eip/T s» 11.8 
at k p TD = 0.022. Analysis of the monotonic solution ip(r) 
at even smaller Debye radii suggests that the magnetic 
potential at the filament center grows approximately as 



eipo/T ~ — A\ak p rD — 3.7, 



(37) 



eip/T = eip /T - 2 ln[(l + (r/r B ) 2 } 
8/(k p r B ) 2 T 



n b 



(1 + (r/rB) 2 ) 2 tnc- 



2 p> 



(39) 
(40) 



where = 2 \/2ru is the radius of the Bennett pinch The 
Bennett distribution is accurate at small distances from 
the filament center, r ~ r^, where magnetic screening is 
negligible. 

In the range of distances re ^ r < k^ 1 , where the 
plasma return current is dominant, one can omit the sec- 
ond term in the right-hand side of Eq . (1521 , 



1 d d / eip 
k 2 r dr dr \ T 



eip 

-jT, 



aK (k p r), (41) 



where a is some constant, and K is the modified Bessel 
function of the second kind. Solutions (f3T)|) and (|4"Tj) 
should match at distances, rs « r « k p l . The mag- 
netic potential in the Bennett pinch has the following 
asymptotics at r >> rs- 



4 In fc„ 



(42) 



when k p rjj — > 0. Therefore, near the filament center (at 
the distances r ~ ro) the plasma return current is much 
smaller than the beam current (etpo/T << l/k 2 r 2 D ). 
On the other hand, because the beam current decreases 
faster than the plasma return current, the latter becomes 
dominant at distances r >> m- At even larger distances 
r > k~ l , these currents reaching the minimum value be- 
come equal to each other. 

This currents' behavior allows us to find the approx- 
imate analytical solution of Eq. (|33|) at small k p ro- 
Omitting the first term in the right-hand side of Eq. 
near the filament center, we obtain 



while the solution (|41[) at r < < k p 1 can be approximated 
by: 



— « a(-7£ + In 2 - lnfc p r), 



(43) 



where je = 0.5772... is the Euler constant. Comparison 
of (|42|) and (|43]l gives a — 4 and the expression for the 
magnetic potential at the filament center, 

^=-41n(VB)+4(- 7B + ln2), (44) 



1 d 

k 2 r dr 



/eip\ 



1 



exp 



>' D 



/eip _ eip \ 
\T T ) 



(38) 



which, after substitution r 2 B — 8r|,, reduces to the ex- 
pression Eq. (I37[) found by the shooting method. Using 
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Eq. (l44l) . we find that the constant n&* in Boltzmann dis- 
tribution at small k p rjj is given by 



nb* 



T exp(-eVJo/T) 



7 2 2 

a x k p r B 



(45) 



where a,\ = exp(47^)/2 w 5.0. At large distances 
from the filament center where the beam and plasma 
currents become equal to each other, the magnetic po- 
tential ipoc is determined by the equation etpoc/T = 
(mc 2 /T)(rib*/n p ) exp (eipoo/T) and equal to 



T 



-W(-aifcprf,) « a\k 2 r 2 Dl 



(46) 



where W(x) is the Lambert W-function. It is convenient 
to use the interpolation formulas for the magnetic po- 
tential and beam particles distribution which work at all 
distances and can be used even for moderately small pa- 
rameter k p rB, 



— = AK a [k p (r + r%) 1 ] + — 



n b 



8/(k p r B ) 2 



T 



(1 + (r/rfl) 2 ) 2 mc 



2 "P 



nbc 



(47) 
(48) 



where ti^oo = (n p T ' /m<r)eil) 00 /T '. These formulas de- 
scribe the modified Bennett distribution of the beam par- 
ticles in the presence of plasma return current. 

We have found a general solution of Eq. (|32|) which de- 
pends on the transverse temperature of beam particles 
T, and the radius of the Bennett pinch tb (or t\d). For 
an arbitrary transverse beam dimension, these two pa- 
rameters can be determined from two conservation laws. 
The procedure is straightforward when the solution of 
Eq. (|32p is found numerically; the result for is presented 
in Fig. [5] by purple bullets. To find analytical formu- 
las, we introduce dimensionless variables tp — eip/T, n — 
mc 2 rib/ '(Trip), and dimensionless parameters T = T/T* 
(where T* = mc 2 nbo/n p ), hoo = mc 2 rn loo / (Tn p ), and 
"00 = e-ipo/T. Substituting Eq. (gSJ) into Eq. (|2"3)l. we can 
evaluate the number of particles as 



8ttT + fi boo Tk 2 L 2 = k 2 L 2 , 



(49) 



The first term in the left-hand side of this equation is 
proportional to the number of the particle in the Ben- 
nett pinch and the second term is proportional to the 
number of particles outside of the pinch. To evaluate the 
potential energy at small k p r B , we use Eqs. l|3"9")) and (159")) 
and transform Eq. (f2"S"]) to 



Tk 2 pL 2 



\f 2 (8^ - 16tt) = -\k 2 v L 2 , (50) 



We can make further simplification, noting that T » 1 
when k p L » 1 and neglecting the term in the right- 
hand side of Eq. (|5"0"|) proportional to the initial potential 




20 40 60 80 100 



FIG. 10: The parameters T and tb of the modified Bennett 
pinch as a function of the transverse beam dimension L. The 
pinch is obtained from Eq. (|32|) by the shooting method (pur- 
ple bullets) and analytically (black lines). The dashed black 
line corresponds to the analytical formula for the radius with- 
out adjustment. 



energy. Using Eq. (1441) and formula hb 
find the temperature and the radius 



T = 



k 2 L 2 



I 



k 2 r 2 - 
K p r B — 



8tt 
ot\k 2 L 2 



. + W(aok 2 L 2 )8^T , 
W[a (k p L + Ai) 2 }. 



otik 2 r B , we 

(51) 
(52) 



Surprisingly despite all approximations, formula (|5"Tj) 
very closely reproduces the result of the numerical in- 
tegration of Eq. ([32]), see Fig. [TO] (left panel). To im- 
prove accuracy of the other formula, we introduce the 
adjustment parameter A\ into Eq. (|52[) . otherwise it is 
not accurate at moderately large k p L. Comparison of 
the adjusted (A\ = 35) and unadjusted (A\ = 0) formu- 
las with the result obtained by the shooting method is 
given on the right panel in this figure. 

Now we find several integral characteristics of the mod- 
ified Bennett pinch. It is interesting to note that the to- 
tal magnetic energy W m in this pinch is a finite quantity 
(in contrast to the regular Bennett pinch, where it di- 
verges at large distances). Using Eq. (|47|) we calculate 
B = — V-0 and then integrate B 2 /8ir over the transverse 
area. The result can be expressed in special functions; 
with some adjustments, it is given at at small k p rs by 



W„ 



T 2 
Aire' 



-(-AAA + 66k p r B - 167rln/c p r B ). (53) 



The number of particles trapped in the modified Ben- 
nett pinch can be found by integrating Eq. (1481) over the 
transverse area. Using Eq. (|51[) and making some adjust- 
ments, we find 



N- 



1 



l + 0.83r B 
+ W(a k 2 L 2 ) ' 



(54) 



where r B should be calculated from Eq. ([52l . Obviously, 
the number of untrapped particles is given by 



N 



W(a k 2 L 2 ) - 0.83r B 



1 



W(a k 2 p L 2 ) 



(55) 
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